'''
import arcgisscripting, os
gp = arcgisscripting.create()

gp.workspace = "f:/redlands"

# Get a list of the feature classes in the workspace.
tbs = gp.ListTables()

# Loop through the list of feature classes.
tbs.reset()
tb = tbs.next()

while tb:
    # Set the outputname for each output to be the same as the input.
    output = out_workspace + os.sep + tb

    # Clip each input feature class in the list.
    gp.AddField(tb, "ECO_CODE", "float")
    tb = tbs.next()
    print "aaa"
'''
'''
from numpy import *
from scipy.cluster.vq import whiten
from scipy import *
x = array([1., 3., 8., 9.])
variance = cov(x, bias=1)                       # normalized by N-1
#print variance
#variance = cov(x, bias=1)
#print variance
mean = mean(x)
y = (x - mean) / variance
#print y
z = whiten(x)
print z

'''
#from scipy.spatial.distance import cdist
import numpy as np

a = np.array([[100,20],[40,50],[70,80],[40,202],[50,20]])
p = 3
a = np.delete(a, p, axis = 0)
b = [[100,20],[40,50],[70,80],[40,202],[50,20]]
del b[p]
print len(b), type(len(a))
'''
XA = np.array([7, 8, 5])
XB = np.array([4, 4, 1])
#c= cdist(XA, XB, metric='euclidean')
dist = np.linalg.norm(XA-XB)
print dist
'''